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Abstract 

We present a semi-decision procedure to tackle first order differential equations, 
with Liouvillian functions in the solution (LFOODEs). As in the case of the Prelle- 
Singer procedure, this method is based on the knowledge of the integrating factor 
structure. 
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1 Introduction 



The problem of solving ordinary differential equations (ODEs) has led, over the 
years, to a wide range of different methods for their solution. Along with the 
many techniques for calculating tricky integrals, these often occupy a large part of 
the mathematics syllabuses of university courses in applied mathematics round the 
world. 

For a first order differential equation (FOODE), finding the solution can be 
equated to determining an integrating factor. A remarkable method for finding such 
factors was developed, in 1983, by Prelle and Singer Q. Their method is based 
on the knowledge of the general structure of the integrating factor for FOODEs 
of the type dy/dx = M{x,y)/N{x,y), with M and polynomials in their argu- 
ments, which present a solution that can be written in terms of elementary func- 
tions (EFOODEs)^. Their approach is very attractive due to the fact that it is 
non-classificatory and of a semi-decision nature. Therefore, it has motivated many 
extensions of the original idea ||2|, ^, | , ^ . 

Using the results presented in |6|, we presented (see 0) a method which is 
an extension to the Prelle-Singer (PS) procedure allowing for the solution of some 
LFOODEsQ {dy/dx = M(x,y)/N{x,y), with M and N polynomials in their argu- 
ments). The method is also based on the general structure of the integrating factor, 
that was concluded to be of the form: R = e"^''^^''^^ YVi^iVi{x,yy\ where rg is a 
rational function of {x,y), the Uj's are irreducible polynomials in {x,y) and the Cj's 
are constants. 

The method presented on 0] used a conjecture about the nature of the Vi^s 
(proved on Q) and was restricted to a class of LFOODEs, namely the ones to 
which ro(x, y) is either f{x) or g{y) or f{x)+g{y), where /, g are rational functions. 
Here, we further detail the structure of the integrating factor thus allowing us to 
remove the above mentioned restriction about tq, maintaining the semi-decision 
nature of the approach. 

The paper is organized as follows: in section ^, we analyze the structure of the 
integrating factor for LFOODEs; in the following section, we show how to apply that 
knowledge to construct a semi-decision method to tackle LFOODEs and show some 
examples of the application of the method. Finally, we present our conclusions. 

2 The Structure of the Integrating Factor 

Based on earlier results [], we will, in this section, improve our knowledge of the 
structure of the integrating factor thus setting the stage for the presentation of a 
semi-decision method to deal with LFOODEs. 
Let us then summarize these initial results: 



^For a formal definition of elementary fmiction, see [^Oj. 

^Liouvillian functions are an extension of elementary functions, see ml 
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2.1 Previous Results 



A seminal result on dealing with LFOODEs was obtained by Prelle and Singer in 
1983 [|. They have demonstrated that, for a LFOODE 

dy ^ M{x,y) 

dx N{x,y) ^ ' 

where M and N are polynomials in (x, y) with coefficients in the complex field C, 
if its solution can be written in terms of elementary functions, then there exists an 
integrating factor of the form = Hi /i^' where fi are irreducible polynomials and 
rii are non-zero rational numbers. Appplying this to FOODEs of the type (||), we 
have: 

i •' * 

where D = Nd^ + Mdy. 

From plus the fact that M and N are polynomials, they concluded that 
D[R]/R is a polynomial and that fi\D[fif\ [|]. We now have a criterion for choosing 
the possible fi (build all the possible divisors of D[fi] up to a certain degree) and, 
if we manage to solve (||) , thereby finding rii , we know the integrating factor for the 
FOODE and the problem is reduced to a quadrature. 

In 1^, 0, next steps were taken: it was shown that, for a LFOODE of type 
(|l|), the integrating factor is of the form: 

n 

R = e-o(^^y)l[v,{x,yr'. (3) 

i=l 

where rg is a rational function of {x, y), the Vi^s are irreducible polynomials in (x, y) 
and the Cj's are constants. 

From this, we could conclude (see [^) that D[rQ] is a polynomial and that the 
Uj's are eigenpolynomials of the D operator. 



2.2 A Theorem Concerning the Structure of rg 

In this section, we are going to demonstrate a result about the structure of rg that 
will allow us to generalize the method presented in 0. In order to do so, we are 
going to use some earlier result^. 

For a LFOODE of the form dy/dx = M{x,y)/N{x,y), where M and N are 
polynomials in {x, y), the integrating factor R is given by i? = e''"^^'^-' IVi=i Vi{x, yY^, 
where ro is rational function of {x,y), Vi are irreducible eigenpolynomials of the D 
operator (where D = Ndx + Mdy),Ci are constants and D[rQ] is a polynomial in 
{x,y)- 



Theorem 1: Let the exponent rg be expressed as P(x,y)/Q{x,y), where P and 
Q are polynomials in {x,y), with no common factor. Then we have that Q\D[Q] 
{i.e., D[Q]/Q is a polynomial in {x,y)). 

Proof: Since D[rQ\ is polynomial (see j^), we can write it as: 

■^In other words, fi is a factor of D[fi]. 
^These results can be found in ||. 
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L»[ro] = D 



QD[P]-PD[Q] 



n 



where 11 is polynomial in Multiplying (Q) by Q, one obtains: 



D[P] 



PD[Q] 

Q 



UQ. 



(4) 



(5) 



Since D is a linear differential operator, with polynomial coefficients, and P is 
polynomial, D[P] is also polynomial. Therefore, since HQ is also polynomial, we 
may conclude that ^^'^^ is polynomial either. Since, by hypothesis, P and Q have 
no common factor, we can infer that (i.e., D[Q]/Q is a polynomial), as we 

wantted to demonstrate. 



This result, in turn, leads to the following corollary: 

Corollary 1: We can write Q as Yli-i Qiix , y)"^\ where the qi's are irreducible 
independent eigenpolynomials of the D operator and the mi 's are positive integers. 



Proof: Q is a polynomial so it can be written as ni=i 9i(^)y)"^S where the g^'s 
are independent irreducible polynomials and the mj's are positive integers. So, we 
can write: 

V j Hi 

implying that gj|-D[gi] (i.e., D[qi\/qi is a polynomial), as we wantted to show. 



3 A Semi-Decision Method to Tackle LFOODEs 

In this section, we will show how to construct, in a semi-decision way (as in the case 
of the PS-method for EFOODEs), a method to deal with LFOODEs of the form 
dy/dx = M{x,y)/N{x,y), where M and N are polynomials in {x,y). 



3.1 Introduction 

We now known that the integrating factor for these LFOODEs is of the form, 

R = eP^^'y)IQ^^'y) S{x,y), (7) 

where we have that: 

• P and Q are polynomials in {x,y). 

• Q can be written as Hi {vqi)^\ where the VqS are independent irreducible 
eigenpolynomials (of the D operator) and the m's are positive integers. 

• S can be written as Wj {'^sj)'^^ ^ where the %'s are independent irreducible 
eigenpolynomials (of the D operator) and the c's are positive integers. 
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In what follows, let us convert this knowledge into a form that will allow us to 
build a solving method. 

Using (|) into D[R]/R = -{d^N + dyM), we get: 



D 

leading to 



+ ^ = -idyM + d,N), (8) 



Remembering that the Vg^s are eigenpolynomials of the D operator, we can write 
-^bsj] — ^sj Vsj, where the As's are the polynomial 'eigenvalues' in associated 
with the %'s. Thus, we can write (P) as: 

i [d[P] - P + Y: As, = - {dyM + d,N) . (10) 

If we then multiply both sides of ( |To| ) by Q, we obtain: 

D[P] - p ^ + Q 5] c, A,,- = -Q {dyM + d^N) . (11) 

3 

Since Q = Oj (^^gi) '"^ ) where the VqS are independent irreducible eigenpolynomials 
of the D operator, we can write: 

where the Ag's are the polynomial 'eigenvalues' associated with the VqS. Using this 
into (^) and re-arranging we get: 

D[P] -PY,mAqi = -\{ {vq^) I E As,' + 9j,M + d^N I . (13) 



This equation will prove to be very important to our method. So, let us now do 
some analysing of its structure: We don't know the degree of the eigenpolynomials Vg 
that build up the function S. Neither do we know the degree of the eigenpolynomials 
Vq that appear on Q. However, once these degrees are set and choosing the values 
for the exponents rrii (thus completing the determination of the degree of Q), we 
can determine the degree of the polynomial P. Let us see how it goes by analysing 
both sides of (|l3|). 

First, note that D[P] = M dyP + N dxP- Since the degree of dxP is at most 
equal to the degree of P minus 1, and the same is valid for the degree of dyP, we 
may conclude that (in what follows, we will denote the degree of any term by dterm): 

^D[P\ < max(dM, djv) + dp - 1. (14) 

Looking now to the second term of the left hand size of ( [I^ ) we may say that 



5 



But we know that D[vqi\ = Xqi Vqi, leading to: 



max(dAg) + d^g < max(dM, d^) + d^g - 1 (16) 

and, consequently, to 

max(d;^) < max(dM, djv) - 1- (17) 
Using (17) into ( |T5| ) we can write 

'^pJ^m.Xq. - max(dM,dAr) +dp- 1. (18) 

Through ( |l4|) and (|l8| ) we see that the maximum degree of the left rand side of (13) 
is max(dM, djv) + dp — 1. Now, looking at the right hand size of (p^), we have that 
its degree is at most equal to dg + max(d;\q, dg^M, dg^jv)- Taking (p^ into account 
and the fact that dg^ < dM — 1 and dg^jv ^ d^r — 1, the maximum degree of the 
right hand size of (|l^) is dg + max(dM, dAr) — 1. Since (|l3|) has to be satisfied, we 
see that the degree of Q bounds that of P. 

3.2 The Steps of the Method 

Based on what has been explained on section (^.l|), let us now stablish the actual 
procedure to be followed. In equation ([l3|) (the cornerstone of our method), we see 
terms involving the eigenpolynomials (and/or the associated 'polynomial eigenval- 
ues') of the D operator (present both on the definition for the function S and on 
for the exponent (Q)). So, the firts job at hand is to determine those. In order to 
do that we have to choose a degree. Let us, for simplicity sake, start with the most 
simple possibility. So, the first step is the 

Step 1 - determination of the eigenpolynomials and of the 'polynomial eigenval- 
ues ' ( of degree 1) of the D operator. 

When these polynomials are determined and put on equation (p^), the next thing 
to do is to 

Step 2 - determine the degree for the Q polynomial. 

Please note that the degree for this polynomial is not fully determined by the deter- 
mination of the degrees for the eigenpolynomials of the D operator, there are still 
the powers for these polynomials (rrij). 

Thus, once the degree for Q is set, we determine all possible values for rrii to 
accomodate this choice with the choice made on the first step just described. So, 
the next step is the 

Step 3 - determination of all possible values for m^. 



After that step is taken, one can see that, on equation (13), we still have P to deter- 
mine. But, from the results present on section ( p.lD , from the previous choices we 
have made, the maximum degree for this polynomial can be determined: we have 
that the maximum degree for the left hand side of (|l3| ) is dp + max(dM, djv) — 1 and 
the maximum degree for the right hand side of ( [T3D is dQ+max(dM, djv) — 1. There- 
fore, once the degree of Q is set, the maximum degree for P is dg + max(dM, djv). 
So we 
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Step 4 - construct a generic polynomial of degree 1 (P = + aix + a2y)- 

Equation ( [l3|) then becomes of the general form: 

P = (19) 

where P is a polynomial. For this to be satisfied, we have to 

Step 5 - equate the coefficients of the different powers of {x, y) to zero. Thus 
generating a set of linear algebraic equations on the undetermined parameters (a's 
and c 's). 

Step 6 - Solve the set of linear algebraic equations thus determining the a 's and 

c's. 

If the solution can not be found, we 

Step 7 - increase the degree of P (up to its maximum degree) and repeat the 
procedure from Step 5. 

If the solution can not be found, we 

Step 8 - increase the degree of Q and repeat the procedure from Step 3 onwards. 

Of course, this procedure can go on indefinitely (remember we do not know the 
bound for the degree of Q). However, the degree of the eigenpolynomials (from 
which we construct Q and S) can also be increased (in Step 1 above it was set to 
be 1). So, instead of trying all the possibilities for Q (they are infinite!), one can 
interrupt this arduous loop and experiment starting all over (from Step 1) but with 
an increased value for the degree of the eigenpolynomials. 

By doing that, we are covering all the possibilities and we may hope that we 
will find a solution within our lifetime. On a brighter tone, most of the examples 
we have come across are solved with low degree for the eigenpolynomials of the D 
operator and for Q. 

3.3 Examples and Results 

In this section, we are going to present examples of application of our method and 
discuss its effectivness. 

First, in order to illustrate the steps of the method just presented, we are going 
to start with a simple LFOODE. This example was artificialy manafuctared by us 
according to two criteria: It is (in a way) simple and it is not solved, as far as we 
know, by any other method. 

Example 1: 

Consider the following LFOODE: 

dy _ {x+l)y 

J ~ 2(2 ^ ' 

ax X — xy — + x'^ 

For this equation, up to degree 1 (Step 1), we have that the eigenpolynomials 
(with the associated eigenvalues) are: 
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• vi = y, Ai = X + 1, 

• V2 = X + y, X2 = 1 + x - y. 

The next step (Step 2) is to choose the degree for the polynomial Q. Starting 
with dg = 1, since we are in the case where d^^ = 1 and d^j, = 1, the only possible 
values for rrii (Step 3) are {mi = l,m2 = 0} and {mi = 0, m2 = 1}. For this 
particular example, the maximum degree for P is 2. So, starting with dp = 1 we 
have P = ai + a2X + a^y (Step 4) and equation ( [l3| ) leads to: 

y™^ {x + y)"'^ {nl {x + I) + n2 {I + x - y) + 3x + 2 - y) + 
(^x — xy — y"^ + a2 + {x + 1) yaS — 
{al + a2x+ aSy) {ml {x + 1) + m2 (1 + x - y)) = 0. (21) 

As we shall see, with the values mi = l,m2 = is possible to find a solution. 
Substituting those into (^) we get: 

y{nl (x + l) + n2 {1 + x - y) + 3x + 2 - y) + (^x - xy - y"^ + x'^'j a2 + 

{x + 1) ya3 -{al +a2x + a3 y) (x + 1) = 0. (22) 

In order to solve the above equation, the coefficients for different powers of (x, y) 
have to be zero. Thus leading to the following set of equations (Step5): 



nl + n2 + 2 = 
-n2 - a2 - 1 = 
-al = 

nl + n2 + 3 - a2 = (23) 
Leading to the solution for the coefficients (Step6): 

al = 0,a3 = a3,nl = 0,n2 = -2,a2 = 1 (24) 
So, the integrating factor for this LFOODE becomes (choosing 03 = 0): 

R = (25) 

{x + y) 

Then the solution is: 

c = y^-l+y)f" _ e-'m^ -^) (26) 

{x + y) ^ y ' 

Our method is designed to deal with LFOODEs, so, in order to analyze its 
effectiveness, we are going to use as our testing ground the LFOODEs found on the 
book by Kamke a traditional testing arena for ODEs. The LFOODEs are the 
equations 1.18, 1.20, 1.27, 1.28, 1.129, 1.133, 1.146, 1.169 and 1.235, as they 
are numerated in the book by Kamke. These equations and their corresponding 
Integrating Factors can be found on table ??. The method deals with all these 
examples successfuly. In order not to make it too cumbrous for the reader, we are 
going to present in detail the calculations for just one of those cases, actualy, the 
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most involved one. This choice for the second example was made aiming to show 
that, even for complex cases, our method is contained and manageable. 

Example 2: 

Consider the following LFOODE (1.169 from the book by Kamke): 

{ax + hf^ + {ax + h){yf + c{yf (27) 

For this equation, up to degree 1 (Step 1), we have that the eigenpolynomials 
(with the associated eigenvalues) are: 

• vi = y, Xi = -yc-by'^ - axy"^, 

• V2 = {ax + b)/a, \2 = ab + a'^x. 

The next step (Step 2) is to choose the degree for the polynomial Q. For this 
particular example, we will see that dg = 4 is needed. Since we are in the case where 
d^j = 1 and djj^ = 1, the possible values for rrii (Step 3) are {mi = 4, m2 = 0}, 
{mi = 3, m2 = 1}, {mi = 2,m2 = 2}, {mi = l,m2 = 3} and {mi = 0, m2 = 4}. 
For this case the maximum degree for P is 8 and the lowest degree to which a 
solution can be found is 4. So, letting P = ai + 02 y + 03 + 04 x + 05 x y + 
aQ x^ + ajx^ + as x^ + ag y^ + aio x^ y + a\\ x y^ + 012 x + ai3 x^ y^ + ou y + ais (Step 
4), equation (13) becomes: 

yTni I j yc — by^ — axy'^j + n2 [ab + a^x j — 2 y (axy + by + c) 

-y^ {ax + 6) + 2 a^x + 2 ab^ + (a^x^ + 2 ax6 + 6^) (2 a2 xy + a4 y^ + a5 y+ 
A a6 x^ + Z a7 x^ + 2 a8 X + Z alO x^y + all y'^ + al2 + 2 al3 xy^) 
-y^ (axy + 6y + c) ^4 ai y^ + x^ + 2 a5 y + 2 xy + a5 x + 3 a9 y^ 
+ alOx^ + 3 all xy^ + 2 al3 x^y + al4^ - i^al y^ + a2 x^y + a3 y"^ + a4 xy^ 

/I O Q Q Q Q 

+ o5 xy + a6 X + a7 x + a8 x + a9 y + alO x y + all xy + al2 x 
+ al3 x^y^ + al4y+ al5^ {ml {-yc - by"^ - axy"^^ + m2 {ab + a^x)) (28) 

As we shall see, with the values mi = 2, m2 = 2 is possible to find a solution. 
Substituting those into (H) we get: 

\Ab^a6 + 6aba7 + 2a^a8 = 0,3b'^ a7 + a'^ al2 + AabaS = 0, 

2aba5 + al2 c+ n2 a^ + 2b'^a2 + 2a^ = 0, -3cai -2ba9 = 0, 
b^al2 = 0,-nl c+ al5b + b^a4 - 2 c = 0,2a'^ a2 + a7c + 6abal0 = 0, 
8aba6 + 3a^a7 = 0,2b'^a8 + 2abal2 = 0,4a^a6 = 0, 
-3b -nl b + b^all - ca3 = 0, a"^ a5 + 3b'^ alO + a8c + Aaba2 = 0, 
bal2 + al5 a + 2b'^al3 + 2aba4 = 0, -nl a-3a-ca4 +2aball = 0, 
ba8 + a^a4 + aal2 + Aabal3 = 0,n2 ab + 2ab + b^a5 + al5c = 0, 
a^all - cal3 = 0, -baS -2ca9 = 0,a6 a = 0,3a^al0 + a6c = 0, 
-2ball - 2aa9 = 0,-3aal =0,-3bal =0,-aal3 = 0,-2aall =0, 

-bal3 - aa4 = 0, -2call - aa3 - ba4 = 0, 

aa8 +ba7 + 2a'^al3 = 0,aa7 + ba6 = 0} (29) 
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Leading to the solution for the coefficients: 



{a6 = 0, all = 0, alO = 0, al = 0, n2 = -1, a2 = 0, a7 = 0, a9 = 0, 

^ ,^c^ -2ah'^al3 ^ „ _ ch 
aS = -1/2 ^ , nl = -3, al4 = ^, 

6^ c 
aiJ = -1/2 — , a5 = — , al2 = -b, a4 = 2 
a a 

al3 = al3, a8 = -1/2 a| (30) 
So, the integrating factor for this LFOODE becomes (choosing as = 0): 

R=^^, rr- (31) 

y3 (ax + b) ^ ^ 



Then the solution is: 



^ - y 2/ (ax + 6) ^^^^ 



2 



4 Conclusion 

In this paper we have present a semi-decision procedure to tackle first order differ- 
ential equations, with Liouvillian functions in the solution (LFOODEs). As in the 

case of the Prelle-Singer procedure, this method is based on the knowledge of the 
integrating factor structure that we now have shown to be: 

n 

R = e'°^'''y^l[viix,yy\ (33) 

i=l 

where vq = P/Q, P and Q are polynomials in {x,y), vi are irreducible eigenpoly- 
nomials of the D operator (where D = Ndx + Mdy), Ci are constants, D[P/Q] is a 
polynomial in {x,y) and (5|-D[Q] (i.e., D[Q]/Q is a polynomial in {x,y)). 
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